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Abstract 

The optimal procedure for extracting the coefficients of different com- 
ponents of a cross section which takes the form of unknown coefficients 
times functions of known kinematical form is developed. When applied to 
e + e~ — * tt+Higgs production at yfs = ITeV and integrated luminosity of 
200 fb _1 , we find that the it — >Higgs CP-even and CP-odd couplings and, 
to a lesser extent, the ZZ -^Higgs (CP-even) coupling can be extracted 
with reasonable errors, assuming the Higgs sector parameter choices yield a 
significant production rate. Indeed, the composition of a mixed-CP Higgs 
eigenstate can be determined with sufficient accuracy that a SM-like CP- 
even Higgs boson can be distinguished from a purely CP-odd Higgs boson 
at a high level of statistical significance, and vice versa. 



1 Introduction 

If Higgs boson(s) (generically denoted as h) exist and are discovered at either 
the CERN LHC or a future next linear e + e~ collider (NLC), it will be extremely 
important to determine both the magnitude and the CP nature of their couplings. 
As reviewed in Ref. |1J, determining the CP properties of a Higgs boson through 
its couplings will be especially challenging. The most promising approaches pro- 
posed to date include: photon polarization asymmetries in 77 — > h B]; momentum 
correlations among the final state r or t decay products appearing in e + e~ — > Zh 
and — > h with h — > t + t~ or it, respectively [0, |[; and weighted cross section 
integrals in pp — > iih at the LHC || and in e + e~ — > ith at the NLC 0. The latter 
ith analyzes examined what can be accomplished using a single observable. Ref. || 
found that under ideal circumstances a SM-like CP-even Higgs boson can be distin- 
guished from a purely CP-odd Higgs boson at a statistically significant level using 
pp — > iih data from the LHC. Ref. found that a statistically significant signal 
for the CP-violating cross term generically present in the e + e~ — > iih cross section 
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for a Higgs boson with both CP-even and CP-odd components might be possi- 
ble. However, neither of these analyzes took full advantage of all the information 
available in the cross section as a function of the kinematical variables. 

In this letter, we outline the optimal technique for determining the coeffi- 
cients Ci appearing in a cross section that can be written in the generic form 
da /dtp = J2i c ifi((p), where (ft denotes the final state phase space configuration. 
The application upon which we shall focus is e + e~ — > tth production, where the 
Q are functions of the Higgs couplings. By extracting the c$ we can determine all 
the Higgs couplings and, thence, its CP nature. This use of the full information 
contained in the final state distributions, as encoded in the c;, leads to significant 
improvement in the statistical precision with which the couplings/CP-nature of a 
Higgs boson can be determined. For example, in e + e" — > tth at yfs = 1 TeV, if 
L = 200 fb _1 and final state reconstruction efficiency is of order e = 0.25, a SM- 
like CP-even Higgs boson can be distinguished from a pure CP-odd Higgs boson 
at roughly the 9.5a statistical level. 



2 General Technique 

We assume that 

where the fi(4>) are known functions of the location in final state phase space, 
0, and the c$ are model-dependent coefficients (taken to be dimensionless in our 
convention). The coefficients Q can be extracted by using appropriate weighting 
functions Wi((f>) such that J = q. In general, different choices for the 

Wi{4>) are possible. However, there is a unique choice such that the statistical error 
in the determination of the c% is minimized in the sense that the entire covariance 
matrix is at a stationary point in terms of varying the functional forms for the 
Wi(4>) while maintaining / Wi((j))fj((j))d(j) = 5ij. Thus, we require 

(a): SVijX j8[w i (<f>)w j (<f>)]J:((f))d(f> = 0, (b): J Sw^fj^) = , (2) 

where Vij is the covariance matrix. The weighting functions which satisfy these 
conditions are of the form 

Wi{</>) = wrj , Wlth x ij = M ij , wh ere M ik = J — — d<f>, (3) 

since, for the Wi(<p) so defined, the constraint (b) implies the minimization condition 
(a) in Eq. @. 

We may then compute c, as 
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It can then be demonstrated that the covariance matrix is 

V l3 = (AqA Cj ) = , (5) 

where <tt = J is the integrated cross section and N = L eS a T is the total 
number of events, with L e ^ being the luminosity times efficiency. The result of 
Eq. (|^) applies only for the optimal weighting functions. 

We note that the above procedure is the optimal one regardless of the relative 
magnitudes of the q. Various limits of the optimal weighting functions have pre- 
viously appeared in the literature. For example, if all of the c$ are small except 
for Cfc, then to isolate the q {i ^ k) the appropriate weighting function reduces to 
Wi{4>) oc fi{<p)/fk{<f)), see e.g. Refs. [|, |, |. 

Our procedure is not altered if cuts are imposed on the kinematical phase 
space over which one integrates. Although such cuts may be required in the actual 
experimental analysis, we have not included cuts in our model computations to 
follow. 

In the e + e~ — > tth process, upon which we shall now focus, in order to fully 
define a point in phase space we must identify the t and t and have no more than 
one invisible particle. Thus, we must employ the final state mode in which one 
t decays leptonically and the other hadronically. Further, the t and t must be 
reconstructed. The overall efficiency for the mixed leptonic-hadronic final state 
decays and the double reconstruction will be denoted by e. Then the effective 
luminosity is given by L e ^ = eL, where L is the total integrated luminosity. We 
shall take e = 0.25. 

If a subset, 0, of the kinematical variables <fi cannot be determined (as is the 
case for e + e~ — > tth if the t,t decay to the purely hadronic or double-leptonic 
final state) the above technique can be applied using the variables, 0, that can be 
observed and the functions fi(<p) = J fi(<f))d<j). This is the case even if one or more 
fk are zero. For example, in e + e~ — > tth the fk{4>) that is a CP-odd function of 
the variables reduces to fa — if one cannot distinguish between the t and t. 



3 Extracting Higgs Couplings in e + e — > tth 

We have applied the above procedure to the extraction of Higgs couplings using 
the process e + e~ — > tth. We define the Higgs couplings via the Feynman rules: 

gm t gmz 

tth: -t{a + zbj 5 )t- , ZZh: c -r^—rg^, (6) 

2m w cos(tV) 

where g is the usual electroweak coupling constant. Thus, a, b, and c are defined 
relative to couplings of SM-magnitude. The SM Higgs boson has a = c = 1 and 
b = 0. A purely CP-odd Higgs boson has a = c = and b ^ 0; the magnitude 
of b depends upon the model — we will display results for 6=1, which would 
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correspond to tan (3 — 1 in a two-Higgs-doublet model of type II (see Refs. |l|, |7| 
for details). This latter choice would, in particular, apply for the CP-odd A of 
the minimal supersymmetric model (with tan/3 = 1). 

The tih cross section contains five distinct terms: E(0) = X)f=i c ifi((p)y where 

C\ = a 2 ; c 2 = b 2 ; c 3 = c 2 ; c 4 = ac ; c 5 = be. (7) 

Of these, the only term in £(</>) that is actually CP- violating is that proportional 
to fee; this is the term upon which Ref. |J focused. Our approach makes use of 
the fact that the full cross section contains additional information regarding both 
b and c. 

We have considered three distinct Higgs coupling cases: 

• I) The Standard Model Higgs boson, with a — c — 1, b — 0. 

• II) A pure CP-odd Higgs boson, with a = c = 0, 6=1. 

• III) A CP- mixed Higgs boson, with a = b = c = 1/ \/2. 

For unpolarized beams, y/s = 1 TeV, m h = 100 GeV and m t = 176 GeV, the 
integrated cross sections in cases I, II and III are o? = 2.71, 0.53, and 1.62 fb, 
respectively. Adopting L eff = 50 fb -1 , we then computed 

X 2 = E ( c i ~ c°)(c, - c°)V , with = *W* , (8) 

i,j=l °T 

(see Eq. (|5|)) as a function of location in a, b, c parameter space, where the c° for a 
given case are computed from the model input values of a, b, c (given above) using 
Eq. (0). Surfaces of constant x 2 = \ and 36 are displayed in Fig. p] for each of the 
three cases. We have indicated the parameter space location of models I, II and III 
by a solid bullet, square, and star, respectively. The x 2 = \ surfaces indicate the 
la errors on the parameter determinations. The \ 2 = 36 (or 6a) surfaces will be 
useful as a reference in assessing the level at which we can distinguish the above 
three model cases from one another. 

Due to the fact that the five q are functions of only the three parameters, a, b, c, 
the x 2 = 1 surfaces in Fig. [I] are not perfect ellipsoids. Nonetheless, we follow the 
usual procedure of defining the ±lcr errors in any one of the a, b, c parameters 
by the largest and smallest values that the given parameter takes as one moves 
about the \ 2 = 1 surface. (These extrema define the locations of the two planes of 
constant parameter value that are tangent to the x 2 = 1 surface.) The resulting 
la errors are tabulated in Table [[]. (The upper and lower limits for a, b, c employed 
for the x 2 = 1 surface plots of Fig. |1] are only just beyond the extrema values.) 
We observe that a is well determined in all cases, but especially for the a ^ 
cases I and III. Similarly, b is well determined in the b ^ cases II and III. The 
magnitude of the error in c is similar for all three cases, and is never especially 
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nput Mode 




P — mixed 




y^=1 surfaces x =^6 surfaces 



Figure 1: Surfaces of constant x 2 = 1 and 36 are displayed for: I) a = c = 1, 
b = 0; II) a = c = 0, b = 1; and III) a = b = c = I/a/2- The parameter 
space locations for I, II and III are indicated by a solid bullet, square, 
and star, respectively. Results are for unpolarized beams, = 1 TeV, 
m h = 100 GeV, m t = 176 GeV and L eS = 50 fb _1 . 
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small. Of course, a much better measurement (e.g. ±5% for a SM-like h) of or 
bound on c will be available from inclusive Zh production; however, this does not 
lead to reduced errors for a and b. Some improvement in the errors is possible if 
the electron beam can be negatively polarized without loss of luminosity; the errors 
for P(e~) = — 1 are given in the table. In what follows, we shall only consider the 
case of unpolarized beams. 

Table 1: We tabulate the la errors, as defined in the text, in a, b, and c for the 
three Higgs coupling cases I, II and III, assuming y/s = 1 TeV, rrih = 100 GeV, 
m t = 176 GeV and L e ff = 50 fb~ x . Results for unpolarized beams and for 100% 
negative e~ polarization are given. 



Case 


Un 
a ± Aa 


Dolarized 
b± Ab 


e 

c± Ac 


F 

a ± Aa 


(O = - 
b± Ab 


1 

c± Ac 


I 


1 +U.043 
-0.066 


r, +0.76 
-0.76 


1 +0.51 

-0.82 


1 +U.U37 

-0.057 


n+0.70 

-0.70 


+U.51 
-0.77 


II 


n +u.iy 

-0.19 


+o.uy3 

-0.14 


n +0.58 

-0.58 


n+U.18 
U -0.18 


+U.U79 
-0.12 


n +U.55 
U -0.55 


III 


1 +U.U75 
-0.087 


1 +0.31 
V2 -0.62 


1 +0.57 
V2 -0.80 


1 +0.U66 
y/2 -0.074 


1 +0.27 
V2 -0.46 


1 +U.56 
-0.73 



Most important is the ability to distinguish different Higgs CP mixtures from 
one another. Referring to Fig. [I], we observe the following:]] 

• If the Higgs is the CP-even SM Higgs boson, then the pure CP-odd case is 
well beyond even the x 2 = 36 surface, and, in fact, it lies on roughly the 
X 2 ~ 90 surface, corresponding to discrimination at the 9.5<r statistical level. 
Even the equal CP mixture case III (the parameter location of which appears 
behind the x 2 = 36 surface in the figure) is ruled out at the 4.8a level. 

• If the Higgs is pure CP-odd, with SM tt coupling magnitude, then the CP- 
mixed and CP-even cases lie 17a and 34<r away, respectively. 

• If the Higgs is an equal mixture of CP-even and CP-odd, with coupling 
strengths specified by a = b = c = l/y/2, then the SM CP -even and pure 
CP-odd cases I and II are both about 6.3a away, i.e. just a bit further away 
than the x 2 = 36 surfaces plotted. 

These results improve if the e~ beam has negative polarization. The discrimination 
abilities are summarized in Table 0. 

x We refer to a parameter location on the x 2 = s 2 surface as an s-sigma deviation in the 
sense that the relative probability or likelihood compared to s = is given by exp[— s 2 /2], just 
as for a one-dimensional parameter space. Thus, % 2 = 36 corresponds to relative probability of 
1.52 x 1CP 8 . This differs from the integrated probability for the parameters to lie outside the 
X 2 = s 2 surface, which for x 2 — 36 is 7.49 x 10~ 8 for 3 parameters, i.e. degrees of freedom. 
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Table 2: We tabulate the number of standard deviations, y/x 1 , at which a given 
input model (I, II or III) can be distinguished from the other two models, assuming 
y/s = ITeV, m h = 100 GeV, m t = 176 GeV and L eS = 50 fb -1 . Results for 
unpolarized beams and for 100% negative e~ polarization are given. 





Unpolarized e 


P(e-) = 


-1 




Trial Model 


Trial Model 


Input Model 


I 


II 


III 


I 


II 


III 


I 




9.5 


4.8 




11 


5.5 


II 


34 




17 


40 




20 


III 


6.3 


6.3 




7.3 


7.3 





We can also analyze our ability to determine that the CP-violating component 
of proportional to C5 = be, is non-zero. We consider model III (the only one 

of our three models for which be 7^ 0). We plot the x 2 = 1 (la) surface in a, b, 
and be space and look for the extrema of be. We find that these extrema occur for 
a ~ b ~ l/y/2 and that be can range from —.05 to +.91, assuming L c g = 50 fb -1 
and unpolarized beams. Clearly, we are not far from establishing a non-zero signal 
at the la level. For twice as much luminosity, L e ^ ~ 100 fb , the extrema of 
be on the la surface are +.15 and +.79, and a non-zero value of be would have 
been established at better than the la level. At the la level, L e fr = 50 fb~ upper 
bounds on \c%\ = \bc\ in models I and II are 0.65 and 0.55, respectively. The above 
results are all somewhat better than obtained for these same models using either 
of the observables (C or (9 ropt ) employed in Ref. ||. 



4 Final Remarks and Conclusions 

In this letter, we have outlined the optimal technique for extracting the coef- 
ficients that appear in a general cross section which is a sum of model-dependent 
coefficients times known kinematical functions. Application of this technique to 
e + e~ — > tih results in good prospects for pinning down the CP nature of the h at a 
1 TeV e + e~ collider operating at an expected luminosity of L = 200 fb -1 , provided 
only that the h has a reasonable production cross section (roughly > 0.5 fb) and 
that the tih final state can be reconstructed with reasonable efficiency (roughly 
e > 0.2). In particular, for a Higgs mass of 100 GeV and unpolarized beams, it 
will be possible to demonstrate, using e + e~ — > tih production only, that a SM 
Higgs boson has the expected CP-even couplings to it and ZZ within and 

-82%' res P ec ti ver y) an d that any CP-odd coupling to ti is less than 76% of the 
CP-even SM strength. The precision with which both the CP-odd and CP-even ti 
Higgs couplings can be determined is somewhat improved for a negatively polar- 
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ized electron beam, assuming there is no loss of luminosity. Most importantly, the 
coefficients of the various terms in the e + e~ — > tth cross section can be determined 
well enough that Higgs CP mixtures that are significantly different from one an- 
other can generally be distinguished at a substantial (sometimes very substantial) 
level of statistical significance. 

We have implicitly assumed that the systematic error in the overall normaliza- 
tion of the tth cross section will be relatively small, e.g. < ±5%. If this is not the 
case, then one can focus on the ratios of the different cross section coefficients to 
one another. Our technique is easily adapted to this situation. 
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